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Abstract: We evaluate the —>■ e~e~^ h process through the ZZ fusion channel at 

the International Linear Collider operating at 500 GeV and 1 TeV center-of-mass energies. 
We perform realistic simulations on the signal process and background processes. With 
judicious kinematic cuts, we hnd that the inclusive cross section can be measured to 2.9% 
after combining the 500 GeV at 500 fb~^ and 1 TeV at 1 ab~^ runs. A multivariate 
log-likelihood analysis further improves the precision of the cross section measurement to 
2.3%. We discuss the overall improvement to model-independent Higgs width and coupling 
determinations and demonstrate the use of different channels in distinguishing new physics 
effects in Higgs physics. Our study demonstrates the importance of the ZZ fusion channel to 
Higgs precision physics, which has often been neglected in the literature. 
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1 Introduction 

The discovery of a 126 GeV Higgs boson at the LHC completes the roster of particles 
predicted by the Standard Model (SM). High-energy experiments now continue their search 
for physics beyond the Standard Model in light of this new era. A major new avenue for 
pursuing this search is the detailed study of the Higgs itself. While the mass of the Higgs 
boson is a free parameter in the SM, its couplings to other particles are dictated by the gauge 
and Yukawa interactions. The observations of this particle are so far consistent with the SM 
expectations, but there is considerable room for new physics to reveal itself in deviations of 
the Higgs properties from the SM. There are also many theoretical scenarios in which such 
deviations would arise at a potentially detectable level. Hence, a precise measurement of 
those couplings is a key tool in establishing a departure from the SM, and in characterizing 
any sign of new physics which may be discovered. 

The LHC will continue to accumulate a large amount of data at unprecedented energies 
for many years, which will improve on the current understanding of Higgs physics. It also 
faces certain limitations intrinsic to a hadron collider, including the uncertainty of large QCD- 
related backgrounds. The LHC can measure particular channels involving specihc modes of 
production and decay in combination, and thus constrain combinations of coupling constants 
and the unknown width. Unfortunately, because it cannot measure a single coupling inde¬ 
pendent of the width, it cannot place strong bounds on the absolute values of couplings, nor 
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on the total width unless additional, model-dependent, assumptions are made [1-6]. Inter¬ 
ference effects can be used to bound the width at a few times its SM value [7-15]. A “Higgs 
factory” such as the International Linear Collider (ILC) has the potential to make precision 
measurements of Higgs physics that take advantage of the simple reconstructable kinemat¬ 
ics and clean experimental environment. One especially appealing feature of the ILC is the 
ability to accurately extract the Higgs width in a model-independent manner. 

The key feature of a lepton collider in making model-independent measurements is the 
ability to determine the inclusive Higgs production rate. This is done using processes such as 
e“e^ h + X where X represents additional measurable particles. Since the initial state, 
including longitudinal momentum, is well known we can infer the Higgs momentum without 
specifying the decay of the Higgs, 


Ph = Pe - e + - Px - (I-I) 

This complete kinematical reconstruction allows us to discriminate the inclusive Higgs signal 
from background and measure the couplings of the relevant production mechanism indepen¬ 
dently of the width. Once this is done, measurements of additional specihc decay channels 
can be used to determine the total width and the absolute values of other couplings. In a 
previous study we discussed this general strategy in detail [16]. Based on available analyses 
the model-independent Higgs width T/^ can be measured at the level of — 5% relative to 
the true width. Most of this error derives from the uncertainty on the inclusive cross section. 
Thus, any substantial improvement of the total width measurement depends critically on 
improving the precision on the inclusive cross section. Currently, the inclusive cross section 
sensitivity is estimated for the “Higgsstrahlung” channel —)■ Zh. The cross section for 

this channel is largest just above the threshold at a center-of-mass energy ^/s ~ 250 GeV, 
where it can be measured using the Z decay to electrons and muons with a relative error 
~ 2.6% [17, 18]. At \/s = 500 GeV the Higgsstrahlung rate is substantially reduced but 
using hadronic decays of the Z may allow one to measure the cross section at ~ 3% [19]. 

Farther improvements can be made by examining the alternate production mechanism 
of ZZ fusion [20, 21] 

e“e''' —>■ e~e^ Z* Z* —)■ e~e'^h, (1-2) 

as depicted in Fig. 1, which has often been neglected in the literature. This mode has a small 
rate at 250 GeV but grows with energy as In^(s/Mf). At 500 GeV it already contributes 
roughly twice as much to the final state e~e~^h as the Higgsstrahlung process Zh —)• e~e^h, 
which falls roughly as I/s, as can be seen in Fig. 2. At I TeV this ratio grows to almost a 
factor of 20. Thus, although the Higgsstrahlung process benefits from a sharp kinematic on- 
shell Z peak through the reconstructible final states into which the Z decays, the ZZ fusion 
channel, which features two energetic forward/backward electrons, should also be exploited 
to make maximal use of the high-energy reach of the ILG. 

In this work we perform a fast detector simulation analysis of the inclusive ZZ fusion 
channel measurement at 500 GeV and I TeV. We simulate the predominant backgrounds 
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Figure 1. Feynman diagram of the ZZ fusion signal process 



4 ^ (GeV) 


Figure 2. Total cross section (in fb) for e e’*' —>■ e e'^ + h at ILC versus ^/s. The dashed curve is 
for Higgsstrahlung mode only. 


and a SM-like Higgs signal and calculate the signal sensitivity using a cut-based analysis and 
multivariate log-likelihood analysis. We find that with the cut-based analysis, we can reach a 
sensitivity on the cross section to the 2.9% level. The multivariate analysis further improves 
the precision of the cross section measurement to 2.3%. 

The rest of the paper is organized as follows: In Sec. 2, we discuss the kinematic features 
for identifying the signal and perform a detailed analyses for the ZZ fusion process at 500 GeV 
and 1 TeV energies including backgrounds. In Sec. 3 we discuss the effects of this additional 
information on the model-independent Higgs width and couplings. We also illustrate the 
potential use of these couplings in constraining higher-dimensional operators. We summarize 
our results in Sec. 4. An appendix is included to address issues relating to potential signal 
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Figure 3. Angles 0i, 02 and (f) as defined in the text. The label e~ (e+) represents the outgoing electron 
(positron) and the Z momentum is given by the difference between outgoing and incoming electrons 
(positrons). The arrows represent momentum directions. The Higgs momentum is perpendicular to 
the plane in the right panel. 


and backgrounds with a single photon in the final state. 

2 Sensitivity Analysis 

We consider the signal process e^e"*" —>■ e~e^h via ZZ fusion as in Eq. (1.2). We assume 
that the incoming leptons are described by the nominal beam energy moving along the beam 
axis in the positive and negative directions respectively. Then the outgoing electrons are each 
characterized by a three-dimensional vector and there are six independent degrees of freedom 
measured in our final state. We choose the dimensionful variables to be the invariant mass of 
the final electron-positron pair rriee and the recoil mass, given by 

T^lec = S - 2yfsEee + (2.1) 

The recoil mass provides the most distinct signal feature since it displays the resonance 
peak at the Higgs mass mh — 126 GeV observable on top of a continuum background. The 
electron-pair mass ruee favors a large value ruge > 250 (600) GeV at a 500 (1000) GeV center- 
of-mass energy. This is distinct from the Higgsstrahlung mode where the pair mass is strongly 
peaked at the Z resonance. Despite a broad distribution for the ee pair mass in the ZZ fusion, 
it still provides some discriminating power against the diffuse electron background. 

The remaining kinematic degrees of freedom can be described by four angles. One of 
these, the azimuthal angle of the Higgs boson around the beam axis, is irrelevant to our 
analysis due to the rotational symmetry of the initial state around the beam line when the 
beam is not transversely polarized. The other three angles, illustrated in Fig. 3, are chosen as 
follows: 6 i is the angle between the intermediate Z coming from the initial electron and the 
Higgs boost direction in the rest frame of the Higgs. O 2 is the angle between the final state 
electron and the Higgs boost direction in the rest frame of the outgoing e“e+ pair. These 
angles take advantage of the scalar nature of the Higgs. The distributions for cos 61 and cos 62 
are rather flat since the Higgs boost direction has no preference to align with the spins of 
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the incoming Zs or ontgoing electrons. There is some correlation between these two angles 
and mild enhancement at larger | cos 0 |, which corresponds to a more collinear configuration. 
This is mitigated by the relatively large virtuality of the Z propagators. In contrast the most 
important backgrounds show much stronger correlation and peaks at high | cos 0 | arising from 
highly collinear regions of phase space which tend to dominate their production. The third 
variable, </>, is defined as the angle between the plane defined by the ZZ pair and the plane 
defined by the outgoing e^e"*" pair when viewed along the Higgs boost direction. It is a 
measure of coplanarity. Here the signal shows a preference for small values of (f), indicating 
coplanar emission of the outgoing e^e"*" pair with the Z propagators and with the incoming 
leptons. This strong correlation is expected since the Higgs does not carry away any spin 
information. The backgrounds will generally have a more complex spin structure which is 
not strongly coplanar. 

In practice, the outgoing electrons of our signal will tend to radiate photons, an effect 
we treat with showering. This radiation degrades our signal resolution. To ameliorate this, 
nearby photons are clustered according to a recombination algorithm and identified with a 
single electron as described in detail in the next section. 

Given our inclusive signal process, the backgrounds are of the form e^e"*" — e~e~^X. 
Obviously, the single photon radiation X = 'j arising from the Bhabha scattering is by far 
the largest. Although the majority of events should be removed by the requirement of a large 
recoil mass mx, beamstrahlung and the effects of the initial-state radiation (ISR), as well 
as the final-state radiation (FSR), will produce additional largely collinear photons. This 
generates a long tail in the recoil mass spectrum due to unobserved photons, mainly along 
the beam pipe. To keep this class of backgrounds under control, we introduce a cut on the 
transverse momentum pT of the outgoing e“e^ pair. Photons which are lost down the beam 
pipe should only contribute small pT differences to the observed final state. Thus the final 
state e“e^ intrinsically has no px as long as collinear photons from final-state showering are 
correctly regrouped with the electrons. The signal, in contrast, has a nonzero px from the 
recoiling Higgs. 

This leaves a background from e“e "*“7 where the extra photon is not close enough to 
either electron to be grouped with it by the clustering algorithm. We find it most convenient 
to simply veto events, in addition to the e“e"*“ pair, with a single isolated photon 

> 10 GeV, > 6 °, ( 2 . 2 ) 

where is the polar angle with respect to the beam. The effectiveness of this cut is illustrated 
in Table 1 for the 500 GeV and 1 TeV runs. (See the next section in Tables 2 and 4 for 
numerical definitions of the cuts.) Simple cuts on invariant mass and px reduce the e“e'’' 7 - 
induced background by 3 orders of magnitude but it remains 30 times larger than our signal. 
However the single photon veto reduces this by more than 90%. 

In principle this affects our inclusiveness. However, the Standard Model processes which 
could produce such a signal, such as /i —)• 77 (where one photon is lost down the beam pipe) 
and h —)• Z 7 , constitute branching fractions of 2.3 x 10“^ and 1.6 x 10“^ respectively. As 


- 5 - 


Cuts (fb) 

Generator level 

^rec? ' f^ee 

PT { ee ) 

Veto isolated single 7 

e~e~^h (500 GeV) 

11.5 

4.11 

3.48 

3.48 

e“e '''7 (500 GeV) 

165000 

317 

67.2 

1.32 

e-e+/i (ITeV) 

24.1 

9.75 

8.49 

8.18 

(1 TeV) 

175000 

1570 

344 

4.73 


Table 1. Cross section (fb) for signal e e~^+h and background e after sequence of cuts. The 
cuts are specified in Table 2 and Table 4 for the 500 GeV and 1 TeV case respectively. 

will be seen, the ultimate precision for the inclusive Higgs production measurement is at the 
~ 2% level so that Higgs decays to 77 or Z 7 would have to be enhanced by more than an 
order of magnitude compared to the Standard Model to be seen in the model-independent 
inclusive measurement. Any such large signal enhancements will be seen at the LHC, to the 
extent that they are not already excluded by current results. See the Appendix for further 
discussion. 

After these cuts some background can remain due to poorly measured final-state particles. 
Particularly at 1 TeV center-of-mass energies, errors on the detected momentum of the final 
state can sometimes fake a recoil mass and a high pT that passes our other cuts. This is nec¬ 
essarily an issue to be determined in detail by experimentalists when working with an actual 
machine and is only parameterized by assumptions on detector smearing and efficiency in our 
simulation. We find that badly measured states are typically associated with very high-energy 
photons. Either these photons are not detected at all due to imperfect calorimeter efficiency, 
or they are reported but with significant error on their transverse momenta. Mismeasured 
low-energy photons will not usually cause a big enough error to satisfy our previous cuts. 
Thus it is useful to veto events with very high-energy detected photons, which are relatively 
rare in the signal. 

Again, one may worry about introducing a bias against photons from Higgs decay, but 
this problem can be addressed. When an event has a high-energy photon we first boost it 
into the rest frame of the Higgs, as determined by the momentum of the outgoing lepton pair. 
If the photon’s energy in the Higgs frame is less than half the Higgs mass, then it potentially 
comes from a Higgs decay, and we do not subject it to the high-energy veto. Thus only events 
with “eligible” photons, 7 *, which could not have come from the Higgs decay, are cut. 

2.1 Simulation framework 

To estimate the expected number of events and derive the sensitivity reach at a given en¬ 
ergy and luminosity we use the ILC whizard setup provided through the detector simulation 
package SGV3 [22]. Beam profiles for several energies have been generated by GuineaPIG [23], 
which includes effects from beamstrahlung and ISR. These profiles are interfaced with whizard 
1.95 [24] to generate parton-level samples. The parton-level samples are then passed to pythia 
which performs showering and hadronization to final-state particles [25] . SGV is a fast detector 
simulation which has been found to agree well with full simulation results. 
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To avoid collinear and soft divergences, at the parton level we require that the energy of 
a final state photon be greater than 10 GeV, and that the invariant masses of final lepton- 
antilepton pairs and of lepton-photon pairs be greater than 4 GeV. We also require that the 
invariant mass of a final-state (anti)electron with an initial (anti)electron, or of a final photon 
with an initial lepton, be greater than 4 GeV. More collinear photons will be generated via 
the showering routines in pythia. 

After simulating tracking and calorimeter hits, SGV attempts to identify charged and 
neutral particles and groups these into jetlike objects according to a sequential recombination 
algorithm. We use the JADE algorithm, which defines a distance between objects 

_ 2EiEj (1 — cos 6ij) ('o o'! 

Vij — p2 ’ 

■^vis 

where Ei and Ej are the energies of two objects and Evis is the total seen energy of the event. 
Nearby objects are merged into subjets until all subjets are separated by yij > 0.01. 

In selecting our observables we first identify the two highest-energy electron/positron 
tracks in an event and discard it if there are fewer than two detected (anti)electrons. We also 
require that these particles have opposite signs. If nearby calorimeter hits included in the 
subjet which contains the track are only identified as photons, then we use the jet momentum 
and energy for our reconstructed lepton. If the subjet contains any particles identified as 
hadrons then we use only the track momentum in order to minimize cases where hadron jets 
overlap with the recoiling electrons. For the purposes of the isolated photon cut described 
above, we define an isolated photon as a jet object which contains only photons and no 
charged tracks or hadronic calorimeter hits. 

In the case of pure photon plus electron/positron backgrounds we simulate both 
e“e +7 and e“e+ —>■ e“e +77 at the matrix element level. After showering there is some 
overlap in the signals described by these two processes. In the spirit of matching calculations 
done for hadron colliders we discard events from e^e"*" —?■ e“e ^7 which produce two isolated 
photons after the clustering procedure. 

2.2 500 GeV analysis 

We proceed with a sensitivity analysis for the ILG running at a 500 GeV center-of-mass 
energy. We apply an initial beam polarization of —0.8 for the electron and -1-0.3 for the 
positron, following the ILG technical design report [26] . We first perform a purely cut-based 
analysis with the cuts listed in Table 2. E* represents only photon hits with energy greater 
than 65 GeV in the rest frame of the Higgs. 

Figure 4 displays the signal and background distributions in rurec, m-ee and the three 
angular variables, after applying Cut 1 as listed. As can be seen, the angular variables show 
considerable distinction from the background which can be used to enhance our sensitivity. 
Cut 2 acts on these angles. 


- 7 - 



Cut 1 

122 GeV < m-rec < 145 GeV 
110 GeV < niee < 370 GeV 
PT{ee) > 40 GeV 
veto 1 isolated photon 

E; < 200 GeV 

Cut 2 

(p < 1.5 


Table 2. Cuts applied at ILC 500 GeV. 







Figure 4. Comparison of signal (solid red) and total background (dashed blue) distributions for 
variables (a) rrirec, (b) niee, (c) Oi, (d) 02 and (e) (j) at ^/s = 500 GeV. Cut 1 in Table 2 is applied. 
For clarity, both signal and background distributions are normalized to unity. 


For this analysis we define the signal sensitivity according to the statistical la relative 










































Process 

Generator level (fb) 

Cut 1 (fb) 

Cut 2 (fb) 

ee —)■ ee/i( Signal) 

11.5 

3.48 

3.11 

ee —)> eeV(,V(, 

659 

23.9 

16.0 

66 G-CU 

78.6 

1.02 

0.70 

ee —)• eeqq 

1850 

9.33 

6.88 

ee —>■ eell 

4420 

5.18 

4.42 

ee —>■ 6677 

1640 

1.18 

0.60 

66 —)> 667 

165 000 

1.32 

0.66 

Total background 

174000 

41.9 

29.2 

5a j a 


8.7% 

8.2% 


Table 3. Cross sections for signal and background processes at ILC 500 GeV. 


Cut 1 

95 GeV < TTirec < 300 GeV 
500 GeV < m-ee < 870 GeV 
PT(ee) > 50 GeV 
veto 1 isolated photon 
< 200 GeV 

Cut 2 

0.14 < 02 < 3.0 
(p < 1.5 


Table 4. Cuts applied at ILC 1 TeV. 


error on the signal, 

5a _ VNs + Nb 
a ~ Ns 

where Ng^b = Lagfi are the expected number of signal and background events after cuts 
respectively. We assume the integrated luminosity L = 500 fb~^ at this energy. The statistical 
significance is then inversely related to the signal sensitivity as Ng/\/Ns + Nb- The effect of 
our cuts on the cross section for signal and background processes is given in Table 3. 

We find that this cut-based analysis can measure the inclusive ZZ fusion signal to a 
relative error of 8%. At this energy the dominant background after our cuts is e~e~^i'ei'e, over 
80% of which is from the process e“e^ —)• W~W~^. The large cross section of e^e"*" —>■ W~W~^ 
is favored by the beam polarization we have used at 500 GeV ILC. It is possible to reduce 
this background with a polarization that favors right-handed electrons; however, this also 
reduces the signal and we do not find any significant gain in sensitivity with the reversed 
polarization. It is possible to enhance sensitivity with an analysis that is sensitive to shape and 
to correlations between variables. This is particularly useful when the signal and background 
display distinct features which are not sharp enough to be efficiently cut on, as in Fig. 4. 


(2.4) 
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Figure 5. Comparison of signal (solid red) and total background (dashed blue) distributions for 
variables Wrec, iTiee, (^ 1 , 02 and (j) at i/s = 1 TeV. Cut 1 in Table 4 is applied. For clarity, both signal 
and background distributions are normalized to unity. 


2.3 1-TeV analysis 

We next extend our analysis to a 1 TeV center-of-mass energy with 1000 fb~^ integrated 
luminosity. The polarization is assumed to be (—0.8, +0.2) as suggested by the Snowmass 
Higgs report [27]. The ZZ fusion process is enhanced with increased center-of-mass energy. 
However, due to radiation from the energetic e~ and e"*", the Higgs mass peak in the nirec dis¬ 
tribution is much more smeared than in the 500 GeV case, and photon radiation backgrounds 
become more significant. The angular variables 02 and cj) show greater distinctions between 
signal and background. To maximize significance we apply cuts as listed in Table 4. 

Figure 5 compares the signal and total background distributions after Cut 1. Table 5 
shows the expected cross sections after Cut 1 and Cut 2. Despite the degradation of the 
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Process 

Generator level (fb) 

Cut l(fb) 

Cut 2(fb) 

ee —)• ee/i(Signal) 

24.1 

8.18 

7.52 

ee —)• eet'e^'e 

978 

31.5 

17.2 

ee eev^^rV^i,,T 

93.9 

3.24 

1.64 

ee —)> eeqq 

2830 

24.1 

13.6 

ee —)> eell 

6690 

13.7 

10.8 

ee —)> eeyy 

3180 

2.68 

1.10 

ee —>■ eey 

175 000 

4.73 

2.28 

Total background 

189 000 

80.0 

46.6 

5a ! a 


3.6% 

3.1% 


Table 5. Cross sections for signal and background processes at ILC 1 TeV with 1000 fb ^ of 
integrated luminosity. 



r 



0.90 0.95 1.00 1.05 1.10 


r 


Figure 6. Five-dimensional Log likelihood as a function of the relative cross section r defined below 
Eq. (2.5) for the 500 GeV case (left) and the 1 TeV case (right). For both analyses, Cut 1 is applied. 


recoil mass peak we gain significance from enhanced statistics and a somewhat improved 
signal-to-background ratio. The cut-based analysis can reach a sensitivity of 3.1%. 

2.4 Multivariate log-likelihood analyses 

To improve upon the cut-based results for reaching the optimal sensitivity, we perform 
a multivariate analysis (MVA) by evaluating a five-dimensional log-likelihood as a function 
of the deviation from the SM. Assuming Poisson statistics in each bin, the log-likelihood is 
defined as 

.^bins 

LL(n; u) = 2 ’^[ui ln{—) + i^i - Ui] (2.5) 

i=i 

where Vi is the expected number of events in bin i for the SM signal plus background, and 
Ui is the number of events in bin i for the SM signal scaled by factor r (signal x r) plus 
background. We evaluate the region around r = 1 and our la deviation from the Standard 
Model value corresponds to ALL = 1. 
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Relative error % 

ILC 250-b500 

ILC 250-b500-bl000 


6.0% 

2.5% 

Improvement 


With HL-LHC 


With HL-LHC 

F 

4.8 ^ 4.7 

4.8 ^ 4.6 

4.5 ^ 3.7 

4.5 ^ 3.7 

9z 

0.99 ^ 0.94 

0.99 ^ 0.94 

0.98 ^ 0.75 

0.98 ^ 0.75 

9w 

1.1 ^ 1.1 

1.1 ^ 1.1 

1.1 ^ 0.89 

1.1 ^ 0.88 

9b 

1.5 ^ 1.5 

1.5 ^ 1.5 

1.3 ^ 1.2 

1.3 ^ 1.1 


Table 6. The improvement on selected coupling precisions by incorporating our ZZ fusion analysis 
from a typical 10-parameter model-independent fit. We show both the ILC exclusive results and ILC 
combined with the optimistic CMS HL-LHC input [27]. For details of htting scheme and combination 
scheme, see Ref. [16]. The results for ILC 250/500/1000 ( GeV) assume 250/500/1000 fb“^ integrated 
luminosities. 

Rather than applying Cut 2 on the angular distributions, we apply Cut 1 and evaluate 
the log-likelihood in the hve dimensional phase space dehned by the variables mreo ^eej 
02, and (/). In the analysis, we perform a 3125-bin analysis by dividing the phase space along 
each variable evenly into hve bins. Figure 6 shows the log likelihood as a function of r. In the 
500 GeV analysis, we hnd the sensitivity on signal cross section improved to 6.0%. For the 
1 TeV case, the multivariate analysis increases the sensitivity to 2.5%. The likelihood prohle 
for the 500 GeV (1 TeV) case is shown in the left (right) panel of Fig. 6. 

3 Impact on Higgs Physics 
3.1 Higgs width and coupling Fits 

Based on our results, the sensitivities on (t],“ which can be reached by studying the ZZ 
fusion channel at 500 GeV and 1 TeV ILC are 6.0% (8.2%) and 2.5% (3.1%) based upon 
MVA (cut-based) analyses, respectively. In combination this yields a 2.3% (2.9%) combined 
uncertainty on from this production mode. 

This is comparable to the current estimated precision of the ILC from studies of Zh asso¬ 
ciate production [19] (that is, cr™'^ of 2.0% achieved by combining 2.6% and 3.0% uncertainties 
from 250 GeV and 500 GeV [28]). Thus, by combining the ZZ fusion and Zh measurements we 
estimate a hnal sensitivity (t 1,“ to 1.5%, a 25% improvement over the Zh channel alone. This 
improvement refines many other derived quantities in the model-independent fit. We demon¬ 
strate the improvement for a few representative quantities in Table 6. We have performed a 
global 10-parameter model-independent ht following Ref. [16]. We compute sensitivities for 
the ILC alone and in combination with projected High Luminosity (HL)-LHC results. We 
take the optimistic projections for HL-LHC precision on cross sections from the CMS detector 
based on Ref. [27]. As discussed in detail in Ref. [16], twice the error of propagates into 
the Ftot determination, and this error dominates for stages beyond the 250 GeV phase of the 
ILG. Our study at the ILG 250-1-500-1-1000 stage relatively improves the total width precision 
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by 16%, Higgs to ZZ coupling by 25%, Higgs to WW coupling by 16%, and Higgs to bb 
coupling by 8%. For other couplings with less precision the is not the largest source of 
uncertainty and less improvement is expected. 

3.2 Operator analysis 

New physics beyond the Standard Model (BSM) could give rise to modifications of the 
Higgs couplings. The proper framework to describe such possibilities in a model-independent 
manner is the effective field theory approach. With respect to the SM gauge symmetry, such 
effects are expressed by dimension-six Higgs operators after integrating out heavy particles 
or loop functions [29-32].^ The operators modifying Higgs to ZZ couplings are naturally of 
particular interest in our case. This is partly because it will be one of the most precisely 
determined quantities through a recoil-mass measurement and partly because it is one of the 
key couplings that could help reveal the underlying dynamics of electroweak symmetry break¬ 
ing. Certain operators may have different momentum dependence and thus measurements of 
differential cross sections may be more sensitive to the new effects.^ The ILC is expected to 
have several operational stages with different center-of-mass energies, and the high-precision 
measurement achievable from ZZ fusion will contribute to our knowledge of these different 
operators.^ 

To demonstrate this important feature, we consider the following two representative 
operators 


Oh = Ohb = (3.1) 

with 

(3.2) 

where cj) is the SM S\J{2)l doublet and A is the new physics scale. The coefficients ch and 
chb are generically of order unity. Following the convention for comparison with existing 
studies [32, 35, 37-40], we adopt the scaled coefficients ch = and chb = -^chb- This 
translates to generic values of ch ~ 0.06 and chb ^ 0.006 for A = 1 TeV. 

The operator Oh modihes the Biggs-ZZ coupling in a momentum-independent way at 
lowest order. This operator renormalizes the Higgs kinetic term and thus modifies the Higgs 
coupling to any particles universally [43, 44]. Equivalently, one may think of rescaling the 
standard model coupling constant. In contrast, the operator Ohb generates a momentum- 
dependent Biggs-ZZ coupling. This leads to a larger variation of the production rate versus 
c.m. energy for the Zh process than the ZZ fusion because of the energy difference in the 
intermediate Z bosons. Consequently, the corresponding deviations of the cross sections are 

'^For recent reviews of these operators, see e.g., Refs. [33-36]. Many of these operators not only contribute 
to Higgs physics, but also modify electroweak precision tests simultaneously [37-40]. 

^For discussions of the effects on Higgs decays due to these operators, see Ref. [41]. 

^Assuming existence of a single operator at a time, limits can be derived, see, e.g., [42]. 
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crb 

Figure 7. Constraints on coefficients of dimension-six operators ch and chb with and without the 
inclusion of the ZZ fusion channel. The dashed and dot-dashed fines represent 2cr deviations from 
zero in the Zh channel at 250 and 500 GeV (blue lines), respectively. The solid (red) fines indicates 
the constraint from ZZ fusion for 500 GeV plus 1 TeV. The outer (black-dashed) contour shows the 
constraint from combined Zh measurements and the middle (yellow) and inner (green) contours show 
the combined 2 a and Itr results with ZZ fusion included. 
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(3.3) 


Such operators receive direct constraints from the LHC from similar production pro¬ 
cesses [37, 38], off-shell Higgs-to-ZZ measurement [45], etc., all of which lack desirable sensi¬ 
tivities due to the challenging hadron collider environment. Based on an analysis of current 
data the coefficient chb is excluded for values outside the window (—0.045,0.075)“^ and ch is 
far less constrained [37, 38]. 


“^The window is (—0.053, 0.044) for single-operator analysis. This smallness of the difference between the 
marginalized analysis and single-operator analysis illustrates that this operator mainly affects Higgs physics 
and thus other electroweak precision observables do not provide much information. 
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We only list above the cross sections which can be precisely measured at different ILC 
stages, with corresponding polarizations taken into account. The distinction between ZZ 
fusion(e“e'''/i) and Z/i-associated production with Z decaying to electron-positron pairs is 
easily made by applying a minimal mge cut above mz- 

In Fig. 7 we plot the expected constraints on the constants Hu and cub from the Zh 
and ZZ processes measured at the ILC, assuming only these two constants among the six¬ 
dimensional terms are nonzero. We show the 95% C.L. contours for different measurements. 
The dashed(dot-dashed) blue line represents the contour from Z/i-associated measurement 
at ILC 250 GeV(500 GeV). The red line represents the contour from combined ZZ fusion 
measurements at ILC 500 GeV and 1 TeV. One can see that at a given energy for a simple 
production mode only a linear combination of the two operators is constrained, resulting 
in a flat direction in the contours. However, measurements of Zh at two different energies 
would allow us to measure both simultaneously, as shown in the gray contour. Moreover, the 
addition of the ZZ information at 1 TeV would offer significant improvements as shown in 
the yellow contour. This allows us to measure cu and cub at the level of 0.04 and 0.004 
respectively. Much of the improvement comes from the fact that va. ZZ fusion, in contrast 
to Z/i-associate production, the Oub operator contributes with the opposite sign of the Ou 
operator. We note here such indirect measurements would strongly constrain BSM physics 
which are otherwise difficult to test, such as singlet-Higgs assisted baryogenesis [46], “neutral 
naturalness” [44, 47, 48], etc. 


4 Conclusions 

To summarize, the ZZ fusion channel for Higgs measurement could provide valuable 
information for precision studies of the Higgs width and couplings because of the logarithmic 
increase of the total cross section versus the center-of-mass energy as seen in Fig. 2. Although 
the signal suffers from large radiation-induced smearing at high energies it can be observed 
with good precision at a 1 TeV run and benefits from a multivariate analysis. We have also 
demonstrated the sensitivity to probe higher-dimensional operators at the ILC, which are 
usually not covered by conventional global fits. We find: 

(i) The inclusive cross section of the ZZ fusion channel can be measured to 2.5% at 1 
TeV. This is competitive with the best estimate of Higgsstrahlung measurement at 250 GeV, 
as shown in Secs. 2.3 and 2.4. 

(ii) Combing the ZZ fusion and Higgsstrahlung channels, the model-independent mea¬ 
surement of the inclusive cross section can be improved to 1.5% with a commensurate im¬ 
provement of the Higgs width determination, as shown in Sec. 3.1. 

(hi) Sensitivities on the inclusive cross section cr^'^ at multiple energies also offer the 
possibility to distinguish contributions from different higher-dimensional operators induced 
by BSM physics. We demonstrate the ability to simultaneously constrain two operators 
whose effects are difficult to observe at the LHC, as shown in Sec. 3.2. Including the ZZ 
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fusion channel provides as large as 50% relative improvement for the constraint on the chosen 
operators compared to the Zh-associated production channel alone. 

In the preceding analysis and discussion, we have shown the appreciable impact of in¬ 
cluding the ZZ fusion channel at the ILC for Higgs physics. Full detector simulations may 
be desirable to further the study of this signal mode. 
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A Consideration of one-photon sensitivity 

As discussed in the main text, we find it useful to simply veto events with a single, 
isolated photon in addition to an electron-positron pair. This cut reduces the potentially 
large background arising from Bhabha scattering plus radiation which can pass the invariant 
mass and pT cuts. This cut also reduces signal events where the Higgs decays to a single 
photon plus invisible particles, or a single photon plus additional particles which are lost 
down the beam pipe. In general we do not expect this to be a relevant effect since our 
final sensitivity for the model-independent cross section is 2.5% while the Standard Model 
processes which might contribute to such events are at the level of 10“^ branching fractions 
or less. Only order of magnitude enhancements to these channels from exotic physics would 
be relevant to our analysis and such enhancements are constrained by exclusive searches at 
the LHG and in future at the ILC. 

Nevertheless, there may be some exotic model which would produce an observable effect 
in the inclusive measurement which is not ruled out by other searches. We note that if one 
wishes to preserve sensitivity to exotic channels which could produce a single isolated photon, 
it is possible to institute cuts which will remove almost all of the background while preserving 
a substantial fraction of any such Higgs decays. We find that, in the reconstructed Higgs rest 
frame, the isolated photon in the background sample is not isotropically distributed. The 
background photon usually appears collinear to the Higgs boost direction, and/or confined to 
be near the radial plane containing the beam and the Higgs boost vector. This is because the 
photon is recoiling against the e“e''' pair with a possible boost along the beam axis due to 
additional unseen photons. We also find that measurement errors on the photon are typically 
larger in the polar angle than in the azimuthal direction. Thus one can largely remove this 
background by cutting on the polar (with respect to the Higgs boost) and azimuthal (measured 
with respect to the Higgs-beam plane) angles of a single extra photon in the Higgs rest frame. 
We find the problematic background can be reduced to the level of a few fb while preserving 
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~ 60% of any hypothetical Higgs decay signal, ^ since the photon from such a decay would be 
isotropically distributed in the Higgs rest frame. Hence any new physics signal large enough 
to affect the inclusive rate would still be observable, although underestimated. 

We note that a cut similar in spirit to this one is already present in the widely used 
analysis of Higgsstrahhmg-inclusive measurement at the 250 GeV ILC [18]. In that case 
additional single photons were removed by a “pr balance” cut when the pT of an isolated 
photon accounted for the bulk of the e“e"*“ pair pT- However, since this more complicated 
approach does not materially change our results we present the simpler case of simply vetoing 
the single isolated photon as described in the main text. 
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